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A versatilidade e o poder de invenção têm 
sido consideradas as características mais 
marcantes da ASEA no tocante a questões 
de electricidade. Por todo o mundo, poderá 
ouvir discussões que terminam com esta 
decisão: "Vejamos o que a ASEA pode fazer!” 
As realizações da ASEA na técnica da energia 
— produção, transmissão e distribuição de 
energia eléctrica — são de renome mundial, 
Pôr a electricidade em uso constitui igual- 
mente uma das esferas de actividade da 
companhia e, aqui também, a pesquisa e 
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“Vejamos o que a ASEA pode fazer” 


desenvolvimento da ASEA têm-se revelado 
duma importância decisiva na criação de 
novos produtos e comodidades, promovendo 
melhorias técnicas nos accionamentos indus- 
triais, automação, fornos eléctricos, tracção 
e transportes. Com facilidades de produção 
em 7 países e representação em 71 países, 
a ASEA forma um grupo de companhias in- 
teiramente internacional em todo o sentido 
do termo. O Sr. achará fácil e útil contactar 
com a organização ASEA para “ver o que a 
ASEA pode fazer”. 


uma força criadora no domínio da electricidade EL ÉCTRICA | LDA. 
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Transmissão de energia 


Tracção e transporte 
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TEMPERATURA, CALOR E ENERGIA 


RESUMO 


Mesno na literatura mais recente é frequente encon- 
trar-se o conceito de calor confundido com o de temperatura 
e ode energia interna, Considerando-se que a estrutura- 
ção cientifica da termodinâmica só se torna possivel desde 
que não se conceda à grandeza calor outros atributos que 
lhe são atribuídos e/ou impostos pelo primeiro princípio da 
termodinâmica, faz-se a análise crítica desses conceitos, 
realçando-se, nomeadamente, a interpretação da célebre 
experiência de Joule e o significado da determinação do 


apelidado equivalene mecânico do calor, 


1) INTRODUÇÃO 


Na grande maioria da literatura relativa à 
Termodinâmica, o calor ocupa indevidamente 
um lugar preponderante, o que se nos afigura 
ser de preferência imposto pela tradição e não 
por razões de ordem científica. Na introdução 
da nossa dissertação (!), salientámos que «a edi- 
ficação científica da termodinâmica só se torna 
possivel desde que haja coragem de romper com 
a tradição, não conferindo à grandeza calor 
outros atributos para além daqueles que lhe são 
reconhecidos e impostos pelo primeiro princípio 
da termodinâmica». 

Afigura-se-nos ser tripla a confusão reinante 
nessa literatura, pela suposta identificação do 
calor com a temperatura, ou com a energia 


(1) Sales Luís — «A Expressão Relativista do Segundo 
Princípio da Termodinâmica» (Tese) — 1965. 
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INTERNA 


por A. SALES LUÍS 
Prof. 1. S. T. 


SYNOPSIS 


Even in wp-to-date literature, the concept of heat is 
usually mistaken for that of temperature and internalenergy, 
But the acientiho framework of thermodynamies is only 
possible when to the heat are not given other properties be- 
sides those related with the firat principle of thermodyna- 
mics. So, a erstical analyses of those concepts is underta- 
ken with an emphasis to the interpretation of the famous 
Joule's experience and the meaning of the determination of 
the currently called mechanical equivalent of heat, 


interna ou com a entropia. A própria terminolo- 
gia, imposta pela tradição, contribui para essa 
confusão; assim, apelida-se calor específico a 
volume constante: 


Cy =[ ado 

oT /Y 
a taxa de variação da energia interna em ordem 
à temperatura, numa transformação considerada 
a volume constante. 

Neste artigo, de naturez+ assás elementar, ten- 
tâmos salientar tais diferenças ao nível do pri- 
meiro princípio da termodinâmica; no artigo 
seguinte, atitude similar é tomada ao nível do 
segundo princípio da termodinâmica. Muito em- 
bora não haja «ideias» novas introduzidas, há 
uma tentativa de «arrumação de ideias», o que 
se nos afigura fundamental para o adequado 
desenvolvimento da termodinâmica. 
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2) TEMPERATURA 


Coloquemos um corpo exposto à incidência 
dos raios solares. Na grande maioria dos livros 
elementares de física, comenta-se o fenómeno 
nos seguintes termos: «esse corpo aquece por- 
que recebe calor do Sol». Em relação a esta afir- 
mação, são pertinentes as duas observações 
seguintes: 


— À experiência revela que a temperatura do 
corpo aumenta; o termo «aquecimento», 
usado na linguagem corrente, significa na 
linguagem científica, aumento da tempera- 
tura tal como o termo «arrefecimento» 
significa diminuição da temperatura. 

— O Sol emite radiações electromagnéticas. 
Como se sabe estas radiações têm um ca- 
ráacter dual: ondulatório e corpuscular. No 
aspecto ondulatório, a energia electromagné- 
tica é transportada pelo apelidado Vector de 
Poynting; no aspecto corpuscular, a energia 
electromagnética é transportada pelos pró- 
prios fotões, partículas constituintes das 
radiações electromagnéticas. Em qualquer 
dos casos, o aumento de temperatura resulta 
da interacção da energia electromagnética, 
que incide sobre esse corpo, com as parti- 
culas constituintes desse corpo. De acordo 
com os actuais conhecimentos científicos, 
não há «calor» transportado pelas radiações 
solares. 


No estado actual da termodinâmica, a tempe- 
ratura é uma propriedade dos sistemas que per- 
mite definir o seu estado térmico, isto é, que 
permite esclarecer se um dado sistema está ou 
não em equilíbrio térmico com outro sistema 
qualquer. De acordo com a sugestão de Fowler, 
pode enunciar-se o apelidado princípio zero da ter- 
modinâmica mos seguintes termos: «Existe uma 
propriedade, apelidada temperatura. A igualdade 
de temperatura é uma condição de equilíbrio tér- 
mico entre dois sistemas ou entre duas porções 
de um dado sistema». 

Embora com um certo carácter arbitrário, a 
temperatura de um sistema, num dado estado de 
equilíbrio, é medida pela indicação instantânea 
de um termómetro. Como se sabe, essa medida 
depende realmente da natureza da substância 
termométrica; no entanto, o segundo princípio 
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da termodinâmica permite definir a temperatura 
de uma forma completamente independente da 
natureza da substância termométrica. 


3) CALÓRICO E CALOR COMO FORMA 
DE ENERGIA 


Na interpretação dos fenómenos térmicos, os 
físicos, no final do século XVIII, consideraram 
pertinente admitir a existência de um fluido, apeli- 
dado calórico, o qual permitia a explicação dos 
fenómenos térmicos então conhecidos. Esses fe- 
nómenos constituem, por assim dizer, os fenó- 
menos inerentes à apelidada calorimetria. Em 
tais domínios, os fenómenos são adequadamente 
explicados em termos desse fluido, calórico ; 
quando se colocam em contacto dois sistemas a 
temperaturas diferentes, o sistema a maior tem- 
peratura cede calórico ao outro sistema até que 
se atinja o equilíbrio térmico. A quantidade de 
calórico cedido por um sistema é igual à quan- 
tidade de calórico recebida pelo outro sistema: 
há conservação do calórico. A diferença entre 
temperatura e calórico pode facilmente ser visua- 
lizada com base na seguinte comparação: a água 
está ao mesmo nível em dois ramos diferentes 
de um sistema de vasos comunicantes, muito 
embora a quantidade da água contida num ramo 
seja diferente da contida no outro ramo; a ccr- 
respondência é a seguinte: 


quantidade de caló- 
“” rico no sistema 


— quantidade de água 
contida no ramo 


— nível —» temperatura 


As conhecidas experiências de Rumford e de 
Joule permitiram evidenciar que, fora do âmbito 
desses fenómenos, o conceito de calórico era 
inadmissível pois conduzia a consequências que 
não estão de acordo com os factos experimentais. 
Saliente-se, no entanto, que — se os fenómenos 
considerados forem exclusivamente os do domí- 
nio da calorimetria — podemos, ainda hoje, inter- 
pretá-los na base da existência do calórico (!). 

No entanto, fora desses domínios, o conceito 
de calórico é inadmissível pois não conduz a uma 


(') Do mesmo modo que os fenómenos electromagné- 
ticos macroscópicos (teoria de Maxwell) podem ser inter- 
pretados na base da existência de um fluido eléctrico. 
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interpretação adequada dos fenómenos térmicos; 
nomeadamente as consequências dela resultantes 
não estão de acordo com os factos experimentais. 

Com a experiência de Joule, aceitou-se que os 
fenómenos térmicos eram adequadamente inter- 
pretados se se admitisse, como induzido dos 
factos experimentais, o princípio seguinte (!), 
dito princípio da equivalência: «Quando, após 
uma série de transformações, um sistema mate- 
rial retorna ao estado inicial, quando, além disso, 
as transformações são tais que o sistema só en- 
tra em interacção com o meio exterior por meio 
de trocas de trabalho e de calor, então: 


a) Se o sistema forneceu calor, recebeu tra- 
balho. 

b) Se o sistema recebeu calor, forneceu tra- 
balho. 

c) Existe uma relação constante entre o traba- 
lho Z das forças exteriores e a quantidade 
de calor posta em jogo. 


Esta relação ] é o equivalente mecânico do 
calor». 

Como se sabe, deste princípio resulta que as 
grandezas trabalho e calor podem exprimir-se 
nas mesmas unidades. Nessas condições, o prin- 
cípio da equivalência traduz-se pela equação: 


t=— ÀQ (1) 


o que, por analogia com a relação entre o 
trabalho das forças conservativas e a diminuição 
de energia potencial: 


= — Amp (1a) 
sugere definir a quantidade (AQ) como a va- 
riação de uma forma de energia que se apelidou 
energia calorífica, energia térmica, ou simples- 
mente calor. Na interpretação dos fenómenos 
térmicos envolvidos em transformações em que 
o sistema retorna ao estado inicial, este conceito 
é perfeitamente compatível com os dados expe- 
rimentais; em tais domínios, pode considerar-se 
o calor como uma energia, (2) como uma proprie- 
dade dos sistemas. 


(+) Jean Becquerel «Thermodynamique», pg. 62, (Li- 
brairie Scientifique J. Hermann — Paris, 1924). 

(*) — Saliente-se, no entanto, que tal interpretação já 
impõe que as forças sejam conservativas. 
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4) GENERALIZAÇÃO: OS SISTEMAS NÃO 
RETORNAM AO ESTADO INICIAL 


Seja, então, um sistema material que se sub- 
mete a uma série de transformações, de acordo 
com as condições impostas no princípio da equi- 
valência, anteriormente enunciado. Considere- 
mos um troço elementar do ciclo descrito e seja 
dz e dQ, respectivamente, o trabalho elementar 
e a quantidade elementar de calor, postos em 
jogo, na interacção do sistema com o exterior, 
nessa transformação elementar: 


dz + dO 


No ciclo descrito, a totalidade de trabalho 
e calor postos em jogo é igual à soma das quan- 
tidades correspondentes a cada um dos troços 
elementares, isto é: 


t+ 3Q = 7 dz + dQ (2) 
Mas, de acordo com a equação (1) é: 
t+ àQ =0 
pelo que: 
f dt+ dQ =0 


Como se sabe, sempre qué um integral cíclico 
é nulo, a integranda é uma diferencial exacta; 
quer dizer, existe uma função U das variáveis 
independentes tal que a sua diferencial total é 
igual a essa expressão diferencial: 


dU = dz + dO (4) 


A generalização do princípio da equivalência 
às transformações tais que o sistema não retorna 
ao estado inicial, impõe a introdução de uma 
nova propriedade, de uma nova função de 
estado U, que se apelida energia interna do sis- 
tema. 

Como se sabe, no caso de sistemas definidos 
por duas variáveis independentes (v, T) tais que 
as únicas forças actuantes se reduzem a forças 
de pressão actuando perpendicularmente à fron- 
teira do sistema, a expressão do trabalho é: 


dtz=—p dV (5) 
e a da quantidade de calor é: 


dQO =Cyv dT+I1IdV (5a) 
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em que Cv e | representam, respectivamente, 
o calor específico a volume constante e o calor 
latente de dilatação. Substituindo (5) e (5a) em 
(4), obtém-se: 


dU =Cv dT+(l—p) dV (6) 


Atendendo a que a condição necessária e sufi- 
ciente, para que uma expressão diferencial seja 
uma diferencial exacta, é que as derivadas cuza- 
das sejam iguais, resulta de (6) a condição: 


Cv: = 9 (|I—p) 
dNV he ( AI a (7) 


o que, atendendo a (5a), evidencia que dQ não 
pode ser uma diferencial exacta. 

Qual o significado físico desta afirmação? No 
aspecto matemático significa que não existe uma 
função Q das variáveis independentes do sistema. Então, 
como qualquer propriedade do sistema tem va- 
lores bem definidos para cada um dos estados 
de equilíbrio do sistema e, consequentemente, é 
uma função das variáveis independentes do sis- 
tema, dizer que dO não é uma diferencial exacta 
significa que O não é uma propriedade do sistema, 
não é uma função de estado. Por outro lado, como a 
energia é, por própria definição, uma função de 
estado, resulta que Q não pode ser uma energia. 

Quer dizer, na interpretação dos fenómenos 
térmicos que ocorrem em transformações não- 
-cíclicas, torna-se inaceitável considerar o calor 
como uma propriedade dos sistemas, como uma 
função de estado, como uma energia. O calor 
é apenas uma grandeza que intervém na interac- 
ção do sistema com o exterior, através da fron- 
teira do sistema. Assim, no caso de um sistema 
definido por 2 variáveis independentes (V, T) 
transitar (*) do estado de equilíbrio AV, T) para o 
estado A + dA (V+ dV, T+ dT), a quantidade 
elementar de calor posta em jogo vale: 


dQ=C,dT+IdV 


Essa quantidade não é definida nem no es- 
tado A nem no estado AL- dA; apenas inter- 
vém quando o sistema evolui de um para o outro 
estado em consequência da interacção do sistema 
com o exterior. Por isso não há variação de 
calor mas calor posto em jogo; por isso, não se 
deve, em relação a uma transformação finita, 


(!) Transformação reversível. 
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escrever 4 Q mas simplesmente Q (calor posto 
em jogo). 

De acordo com o estado actual da termodinã- 
mica, considera-se que os factos experimentais 
são adequadamente interpretados, admitindo-se o 
seguinte princípio, apelidado primeiro princípio 
da termodinâmica: «Existe uma função de es- 
tado, apelidada energia interna. Numa transfor- 
mação elementar, em que a interacção com o 
exterior se reduz a trocas de trabalho e de 
calor, através da fronteira do sistema, a variação 
elementar da energia interna é igual à soma algé- 
brica do trabalho elementar e do calor elementar 
postos em jogo nessa transformação elementar : 


dU=dz+dQO» (8) 


Como se sabe, qualquer sistema material é 
constituído por átomos e moléculas os quais 
estão em movimento e entre os quais se exer- 
cem forças de interacção. Deste modo, o sistema 
material possui energia potencial e energia ciné- 
tica (microscópicas); a sua totalidade corres- 
ponde, a menos de uma constante adicional (!), 
à energia interna do sistema. 


5) INTERPRETAÇÃO DA EXPERIÊNCIA DE 
JOULE 


Nos livros elementares de física, a célebre ex- 
periência de Joule é descrita e comentada em 
termos que a seguir tentaremos reproduzir (*). 
O aparelho é, por exemplo, constituído por um 
conjunto de pás, montadas sobre um eixo; esse 
eixo é posto em rotação por um conjunto de pe- 
sos e, desse modo, as pás são forçadas a rodar 
imersas num líquido (por exemplo, água). Co- 
menta-se: «(Os pesos realizam trabalho mecã- 
nico e este trabalho é convertido em calor pelo 
atrito no líquido. Como o líquido está imerso 
num calorímetro, a quantidade de calor gerada 
pode ser calculada a partir do aumento resul- 
tante da temperatura de tal modo que pode ser 
determinado o equivalente calorífico do trabalho 
mecânico realizado pelos pesos». 

Saliente-se a afirmação: «o trabalho é conver- 


(') Esta constante adicional só intervém se houver 
transformação de matéria em energia e vice-versa, o que 
é estranho à termodinâmica clássica. 

(?) Ver, por exemplo, R. Kronig, «Texthook of Phi- 
sics», pg. 63 | Pergamon Press). 
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tido em calor» e repare-se que não há qualquer 
fundamento experimental para essa afirmação. 
As pás rodam no líquido e a experiência revela 
apenas que há aumento da temperatura do lí- 
quido. Como veremos no artigo seguinte, esse 
aumento de temperatura resulta da produção de 
entropia; no entanto, esta afirmação não é direc- 
tamente induzida da experiência mas é conse- 
quência do segundo princípio da termodinâmica. 
De acordo com o primeiro princípio da termo- 
dinâmica (equação 8): 


dU=dz+dO (8) 


resulta que, estando o calorímetro isolado do ex- 
terior, é: 


dQ=0 
e, portanto: 
dU = dz 
ou 
AU = (9) 


Quer dizer, de acordo com o primeiro princi- 
pio da termodinâmica, o trabalho mecânico rea- 
lizado pelos pesos é fornecido ao líquido e pro- 
voca aumento da sua energia interna. Esse 
aumento da energia interna é acompanhado pelo 
aumento da temperatura do líquido. Os comen- 
tários apresentados tornam evidente a confusão 
que reina entre calor e energia interna; na ver- 
dade, o trabalho não é convertido em calor mas 
provoca o aumento da energia interna do sis- 
tema. Saliente-se, por outro lado, que a expe- 
riência descrita não está nas condições impostas 
pelo princípio da equivalência pois o sistema 
não retorna ao estado inicial. 

Consideremos, no entanto, uma experiência 
complementar da anteriormente descrita, em que, 
utilizando o mesmo calorímetro, se provoca a 
passagem do líquido do estado inicial i ao estado 
final, ambos de equilíbrio e idênticos aos consi- 
derados na experiência anterior, mas de modo 
que haja apenas «calor» fornecido ao sistema 
(z == 0). Neste caso, a equação (8) escreve-se: 

AU =Q (10) 

Então, como a variação de energia interna é a 
mesma nas duas experiências realizadas, resulta 
de (9) e (10) que: 


t=Q (11) 
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Saliente-se, portanto, que o trabalho não é 
convertido em calor. A experiência revela apenas 
que é equivalente o trabalho Z e o calor Q, 
postos em jogo em duas experiências diferentes, 
na medida em que, separadamente, conduzem 
ao mesmo valor da variação de energia interna ; 
numa experiência só há trabalho realizado e na 
outra só há calor fornecido. 

Como salienta Planck, (!) o primeiro princípio 
não especifica apenas que a variação da energia 
interna é igual à soma algébrica do trabalho e 
do calor postos em jogo na transformação ; espe- 
cifica também que essa soma algébrica (isto é, 
a variação da energia interna) é independente 
do modo como se processaram as transforma- 
ções que conduzem o sistema do estado inicial 
ao estado final. É nesse sentido, que se pode fa- 
lar no equivalente mecânico das acções exterio- 
res. Por outro lado, a elevação da temperatura 
de um líquido pode efectuar-se de 2 modos dis- 
tintos: ou fornecendo-lhe directamente calor ou 
por atrito, como na célebre experiência de Joule. 
E esse facto que permitiu a Joule calcular o equi- 
valente mecânico do calor. 


6) CONSIDERAÇÕES FINAIS 


Muitos outros exemplos podiam ser citados no 
sentido de evidenciar que, ainda hoje, reina 
grande confusão entre calor e temperatura e en- 
tre calor e energia interna. É frequente encon- 
trar, na literatura relativa à teoria cinética dos 
gases, que o calor corresponde à energia ciné- 
tica das moléculas. Há, nesta frase, uma dupla 
confusão : por um lado, é a temperatura que for- 
nece uma medida indirecta da energia cinética 
média das moléculas; por outro lado, essa ener- 
gia cinética é uma parte da energia interna do 
sistema. 

As considerações anteriores são relativamente 
elementares e resultam naturalmente do pri- 
meiro princípio da termodinâmica. No entanto, 
afiguram-se-nos fundamentais para uma ade- 
quada estruturação científica da termodinâmica 
e desempenharam papel preponderante na for- 
mulação relativista do segundo princípio pro- 
posta na nossa dissertação. 


(!') Max Planck, «Leçons de Thermodynamique» (tra- 
dução, Librairie Scientifique A. Hermann — Paris, 1913), 
pg. 292. 
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|| SEMANA DE ENGENHARIA QUÍMICO-INDUSTRIAL 


Está em adiantada preparação a II Semana de Engenharia Químico-Industrial que, como noti- 
ciâmos, se realizará no Instituto Superior Técnico, entre 10 e 21 de Abril próximo. 
A exposição estará dividida em quatro partes: 


— Parte industrial — que terá lugar no átrio do I.S.T. 

— Parte documental — contando com a colaboração de vários organismos oficiais e a realizar 
na galeria sobre o átrio do 1.5.T. 

— Parte bibliográfica — com a participação de várias livrarias técnicas e dos grupos de estudo 
de algumas empresas 

— Equipamento de laboratório — a realizar no salão nobre. 


Além de sessões de filmes técnicos e de visitas de estudo, está planeada uma série de confe- 
rências sobre vários assuntos relacionados com a Engenharia Química. Destacamos: 


Dia 10, às 10 horas — «Relações entre a Escola e a Indústria», pelo Sr. Eng.º José Mercier 
Marques, Bastonário da Ordem dos Engenheiros. 


Dia 11, às 10 horas — «A investigação e o desenvolvimento da Indústria Química em Portu- 
gal», pelo Sr. Prof. Dr. Alberto Ralha, Director do Laboratório da Poli- 
cia Científica. 


Dia 14, às 10 horas — «Problemas de Corrosão — Protecção de Materiais», pelo Sr. Eng.º 
Sidónio Geada, Director do Núcleo de Electroquímica e Corrosão do 
IN.LI. 


Dia 15, às 10 horas — «Aplicação de Computadores na Indústria Química», por um técnico 
francês da I.B.M. 


Dia 17, às 15 horas — «Termodinâmica Química», pelo Sr. Prof. Eng.º Sales Luís. 


Dia 21, às 18 horas — «Panorama actual da Indústria Química Portuguesa», pelo Sr. Eng.º 
Carlos Alves, Presidente da A,I.P. 


Realizar-se-ão ainda mais três palestras, cuja efectivação carece ainda de confirmação. 

Nos dias 12 e 19, à tarde, realizar-se-ão dois colóquios que versarão os seguintes temas: «Es- 
truturação do Curso de Engenharia Químico-Industrial» e «O problema dos estágios no Curso de 
Engenharia Químico-Industrial». 

No dia 21, após a sessão de encerramento, realizar-se-á um jantar de confraternização entre 
antigos e actuais alunos, e professores e assistentes do I.S.T. 

O programa definitivo será em breve tornado público. 
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PONTES ROLANTES, GUINDASTES E 


"CONDENSADOR DE SUPERFÍCIE (7.300 m) PARA ÁGUA SALGADA, DO 


TURBOGR 


UPO | DA CENTRAL DE CARREGADO. 


PESO TOTAL (c tubos de refrigeração ) = cerco de 180 + 
COMPRIMENTO = cerca de 14.2 metros 


LARGURA 
ALTURA 


= cerco de 5.5 metros 
= cerco de 7.2 metros 
- (fabricado sob licenço BROWN BOVERI CIE)* 


APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico . 
Fabrico segundo licença de A.C. M. de Vevey, S. 4. 
Fabrico segundo licença de Brown Boveri, Cie. 

Projecto e fabrico segundo licença de Foster 


TURBINAS HIDRÁULICAS — | 
TURBINAS A VAPOR 
CALDEIRAS A VAPOR 


EQUIPAMENTOS E INSTALAÇÕES 
INDUSTRIAIS 


Wheeler, Co. 


CONSTRUÇÕES METALOMECANICAS MAGQUE :..-. 


ALVERCÇCA DO RIBATEJO - PORTUGAL 


TECNICA XII 


O a dna 


Para 0s motores mais aperfeiçoados, o mais aperfeiçoado dos lubrificantes. 


PERRINO APROVA BP VANS 


Os novos motores Perkins foram projectados 
para o dia de hoje; — com [lexibilidade; para 
satisfazer a dupla exigência de alias velocida- 
des nas auto-estradas e cóniinuas paragens e 
arranques no intenso trânsito citadino. Em 
quaisquer condições eles aliam à eficiência à 
economia. Como o mesmo sucede com BP 
Vanellus, a Perkins aprovou-o para todos os 
seus motores. BP Vanellus conserva Os moto- 


res mais limpos, livres de depásitos nocivos 
e de “sludge” Permite aumentar os perio- 
dos entre mudanças de Óleo. A imedials 
eprovacão do BP Vanellus pela Perkins é 
singl evidente do seu maior avanço como 
lubrificantes de veiculos comerciais Os 
proprietários de frotas, para maior rentabili- 
dade do negócio, precisam de BP Vanellus 
para os seus motores (diesel e gasolina). 


O ÓLEO MAIS AVANÇADO PARA FROTAS COMERCIAIS 


TECNICA XIV 


VANELLUS 


C. D. U. 624.374.3 


CIRCUITOS LÓGICOS DE RESISTÊNCIAS — TRANSISTOR 
DIMENSIONAMENTO * 


SUMÁRIO 


No presente artigo expõe-se um método fácil de dimen- 
sionamento de circuitos lógicos de resistências-transistor. 
Tomam-se em consideração as tolerâncias dos componentes 
e das tensões de alimentação. Definem-se e determinam-se 
os factores de entrada e saída necessários para posteriores 
projectos de sistemas lógicos, 


1 — CIRCUITOS LÓGICOS DE RESISTÊNCIAS 
TRANSISTOR 


Salvo o inconveniente do tempo de propaga- 
ção, que limita a sua aplicação a sistemas lentos, 
o circuito lógico de resistências-transistor con- 
duz em geral a grande economia de componen- 
tes, em relação a outras soluções. Não será neste 
artigo, analisado o regime transitório do sistema, 
mas considerar-se-ão apenas as condições míni- 
mas de garantia dos estados de saturação e corte. 

O circuito básico (fig. 1) forma um bloco ló- 
gico de M entradas e uma saída. 

O transistor funciona na região de saturação, 
na região de corte, ou transita muito rápida- 
mente de uma para outra através da região activa, 
Em geral é desprezável a potência média dis- 
sipada na junção colector-base, pois nesta junção 
a corrente ou a tensão têm valores muito baixos, 
podendo tomar-se para temperatura de junção 
o valor da temperatura ambiente. A energia é 
então praticamente toda dissipada nas resistên- 
cias Rc, Rm e Rym. do circuito. 

Pressupõe-se que cada uma das M entradas 
se encontram ligadas a saídas de circuitos seme- 
lhantes, de modo que o estado de cada transis- 
tor fique determinado pela combinação dos esta- 
dos daqueles que imediatamente o antecedem. 
Assim, o estado de saturação dum transistor ve- 
rifica-se quando algum dos transistores ligados 


por LUÍS VICENTE 


Engenheiro Especialista 
do Laboratório Nacional 
de Engenharia Civil 


SYNOPSIS 


In this paper, a straightforward design procedure oy 
RTL circuits is explained, Resistor, transistor and bias 
voltage tolerances are taken into account. Fan in and fan 
out factors are computed. 


à entrada se encontra no estado de corte. Por 
outro lado o corte dum transistor resulta do facto 
de todos os transistores de entrada se encontra- 
rem saturados. 


Vec 


M ENTRADAS 


YBB 
Fig. 1 — Circuito básico 


Associando a estes estados 2 valores lógicos, 
por exemplo : 


Corte — 1 
Saturação — O 
tem-se: 


S=E+E +... .+ Ey 


(*) Estes circuitos pertencem à classe NEM. As suas possibilidades sob o ponto de vista lógico foram descritas, 


duma forma geral, na «Técnica» n.º 347, páginas 479 a 482. 
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em que S representa o estado de saida e Ei, 
Es», ... Em os estados de entrada (!). 
Pela 2.º lei de De Morgan obtém-se ainda: 


S = E. Es... Em 


isto é, a conjunção das negações. É nesta rela- 

ção que se fundamenta a designação de blocos 

NEM atribuída a estes circuitos (Fig. 2). 
Poder-se-ia ter convencionado : 


Corte — O 
Saturação — 1 
pelo que viria: 
S=Ej Es... Em ou 
S=E+E +... +E (1.º lei de De Morgan) 


Fig. 2 — Bloco NEM 


No entanto a 1.º convenção é a correntemente 
utilizada. 


2 — PARÂMETROS, TOLERÂNCIAS 
E SIMBOLOGIA 


No que respeita a transistores de comutação, 
os fabricantes, em geral, especificam o ponto de 
funcionamento no estado de saturação, definido 
pelas correntes de colector Ic e de base Ip. Para 
este ponto são dados: 


Vere — Tensão base-emissor máxima 
Vcem — Tensão colector-emissor máxima 
Vcen — Tensão colector-emissor mínima 


No caso do valor do último parâmetro não 
ser fornecido convirá considerar-se: 


Vcen = 0 
Estes valores são característicos de cada tipo 


de transistor e constituem limites de dispersão 
(Fig. 3 e 4). 


(') O traço representa negação e o sinal - disjunção 
inclusiva (V. «Técnica» n.º 347, pág. 479-482). 
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Quanto à corrente de fuga colector-base, cos- 
tuma fornecer-se a lei de variação (exponencial) 
com a temperatura de junção e o seu valor má- 
ximo para uma dada temperatura (em geral 25ºC). 
Fixado o valor da temperatura ambiente máxima, 
é fácil determinar o máximo valor da corrente 
de fuga, que designaremos para Icou. 


YCE 


YCEN YCEM 


Fig. 3 — Tensão colector emissor 


Fig. 4 — Tensão base-emissor 


Cada equação que a seguir vamos estabelecer 
corresponde a uma condição necessária de fun- 
cionamento. Nela entrarão os parâmetros com os 
valores mais desfavoráveis a essa condição. 

No caso das resistências, há que considerar a 
sua tolerância. Se Rn é o valor nominal duma 
resistência, o valor efectivo Re estará compreen- 
dido no intervalo: 


Ry (1-7), Rx (1 +7) 


TÉCNICA N. 364 


Se o valor nominal duma resistência é fixado 
«à priori», entrará na equação o limite deste 
intervalo mais desfavorável à condição de fun- 
cionamento : 


Rs (1—7) ou Rx (1+7) 


Se R é um valor calculado, o valor nominal 
escolhido Rw pertencerá ao intervalo: 


pois o espaçamento dos valores nominais em 
cada série de resistências é aproximadamente 
igual ao dobro da tolerância dessa série. 

Como o valor efectivo Re existe, em relação 
ao valor nominal Rvw, no intervalo: 


Ry (1-7), Rs (1+7) 


ele existirá, em relação ao valor calculado R, no 
intervalo: 


RA-RA+? 
ou, atendendo a que em geral y << 1: 


R(1—-27,R(1+27) 


Apresentamos a seguir uma lista dos símbolos 
adoptados, com indicação do sinal das correntes 
e das tensões correspondentes ao caso NPN. 
Quando se utilizam transistores PNP todos os 
valores mudam de sinal. 


Is — (-+)-— Corrente de base-saturação 

lc (+) — Corrente de colector-saturação 

Icom (+) — Corrente máxima de fuga de co- 
lector-base 

Iu (+) — Corrente mínima por entrada-sa- 
turação 

lo (+) — Corrente mínima de saída-bloqueio 

VB (—) — Tensão de alimentação 

Vegm (—) — Valor máximo de Vez | Vou | 


| < | Ves| 
Venn (—) — Valor mínimo de Vps. | Vs | 
| <<] Vesl 


Ves (+) — Tensão emissor base-bloqueio 

Vpe (+) — Tensão máxima base-emissor-sa- 
turação 

Vcc (+) — Tensão de alimentação 

Vccm (-+) — Valor máximo de Vcc. | Vccm | 


| > | Vec| 
Vcen (+) — Valor mínimo de Vcc. | Vcen | 
| < | Voecl 


Vce (+) — Tensão colector-emissor-bloqueio 


TECNICA N.º 364 


Vcem ( +) — Tensão máxima colector-emissor- 


-Saturação 
Vcen (+) — Tensão mínima colector-emissor- 
-Saturação 
Re Resistência de carga de colector 
Rm — Resistência de entrada 
RrTM — Resistência de polarização de base 
Y — Tolerância das resistências 
M — Número de entradas 
Em — Factor de entrada 
Fo — Factor de saída 


A, B, I (+) — Parâmetros auxiliares. 


3 — CONDIÇÕES DE FUNCIONAMENTO 


a) Condição de saturação 


Yec 


YcEN 


a === 


YCEN 


Fig. 5 — Condição de saturação 


A condição de corte de um só transistor 
dos M antecedentes deve ser suficiente para im- 
por ao transistor em causa a condição de satu- 
ração, mesmo que os outros M—1 transistores 
se encontrem saturados. A corrente mínima sufi- 
ciente para a saturação é dada por: 


Vez — VBBm 


à — VBBm |, VBE— VCEN (M 
1—27%) Rrm 


—1) 
1—27) Rm 


Im=Ip + 


O valor de Rm deve garantir |m através da 
equação : 


a Vce — VpE 
PR 
(1+27)) Im 


(1) 


b) Condição de corte 


Um transistor deve cortar quando todos os 
antecedentes saturam. Nestas condições a tensão 
emissor-base deve exceder um dado valor Vep: 


CEM 
YcEM 
Meia = === use 
YcEM 
YBB 
Fig. 6 — Condição de corte 
V 
Vcem + VeB ,y Ligas — Ves + VeBN 
(1—27%) Rm (1 + 27) Rrm 
donde se tira: 
 U 
Ves + Venn 


Rr = = vw 
Ve - V 
SEM TE M$ Icom (+27) 
(1—27) Rm 


Vep arbitra-se em geral com um valor posi- 
tivo (transistores NPN): 20º/, de Vpp, por exem- 
plo. O aumento de Vep contribue para uma re- 
dução do tempo de armazenamento, embora com 
sacrifício do valor de Max. 


c) Factor de entrada 


Eliminando Rm e Rrm nas equações anterio- 
res, obtém-se: 


(3) 


sendo: 


a (Ver — Vbsn) q +24) 


(4) 
(Ves + Vsn) (1 — 27) 


Icom 
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fu ja E (Ve — Vcen) (1 + 27) 
(Vce + Vse) (1—2 7) 


(5) 


(6) 


o mm 


EP Nic ) (Vem E Vin +25) 
(Veg ++ Vess) (Ver —Vpe)(1—27)? 


B=A-l— 


e m— 


Verifica-se que Im é tanto menor quanto me- 
nor for M. Para M =1, tem-se 


Iu = lh = a 


A-B 


Factor de entrada para um dado número de 
entradas M é definido pela relação : 


ly A— B 


Ey =— = ———— 
lh A — MB 


(7) 


Para M = 1, vem, necessáriamente, Fi = 1 


d) Factor de saída 


Se fixarmos o valor de Vce, a corrente mínima 
de saída lo e, nessas condições : 


— Vcen— Vek 


Re (1 + 7) 


O 


— Icom 


Yec 


Fig. 7 — Transistor cortado 
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O novo Nível automático 
de engenheiro 


WILD NAK 2 


com circulo horizontal 


ja fimo, 


HEERBRUGG 


Para nivelamentos de média e grande precisão; 
com micrómetro de lâmina de faces plano- 
-paralelas para nivelamentos de grande precisão. 


Erro médio para 1 Km de 
nivelamento duplo: + 1,5 mm 


Erro médio para 1 Km de 
nivelamento duplo com micró- 
metro: + 0,4 mm 


Luneta de 30 x de aumento 
com imagem direita de obser- 
vação 


Solicitem catálogo ou demonstração 
aos representantes exclusivos 


WILD PORTUGAL, LDA. oi ttieon ano 
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ESTAMPARIAS TINTURARIAS 
LAVANDARIAS 


PARA OS VOSSOS TECIDOS TÉCNICOS: 
TRANSPORTADORES, MANCHONS, LAPPING, WOLFRIES, 
FLANELAS, CANEVAS, ETC,, 

EM FIBRAS NATURAIS OU SINTÉTICAS 
CONSULTEM: 


FANAF EL 
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Sendo Rc um valor normalizado que se de- 
termina por: 


Ro = “CC 


lc 


Fácilmente se verifica que a potência máxima 
é fornecida à carga quando: 


Vccen— Re (1 + 7) Icom 
2 


Vce = 


e como, em geral, Rc (1-7) Icom << Vccn, fi- 
xa-se: 


Vcen 
Veg = —— 


: (9) 
Então: 

Ei Vce 

ii a 

Rc (1 +7) 

O factor de saída será, por definição: 

lo (A-B) Vce 

Pisa EO q ms RES 
h  Re(14)1 (80) 


e) Acoplamento de blocos 


No acoplamento de blocos deve observar-se 
sempre a condição: 


Fo > & ny Fm 


sendo ny o número total de entradas de circuitos 
de M entradas. Por outras palavras: a soma dos 
factores de entrada relativos a todas as entradas 
que ligam a uma saída deve ser igual ou menor 
que o factor de saída (Fig. 8): 


Fo > R+2 E + E 


Por outro lado vimos que o factor de entrada 
aumenta com M: 
E A—B 
A-M B 


O seu valor máximo não pode exceder Fo, 
pois na hipótese mais favorável duma saída ali- 
mentar uma única entrada, se tem: 


Fo >, Fm 
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ou seja 
XA—B) Ver ., A-B 
Re (1+7)1 A — MB 
donde se conclui que: 
Mmax <. = eis ai (11) 
N Fo j 
sau" Ni 
/ 
F3 
Nº 
/ 
é 
Nd 
/ 
F, N ' 
/ 


Fig. 8 — Ligação de blocos 


4 — DIMENSIONAMENTO. EXEMPLO 


O número de circuitos necessários para a reali- 
zação de um dado sistema lógico diminui quando 
Fo aumenta. Facilmente se verifica que, para um 
dado tipo de transistor, Fo é tanto maior quanto 


I 
maiores forem | Vcc|, |Ves | e a - Caso não 
B 


haja outros condicionalismos devem sempre uti- 
lizar-se as tensões de alimentação máximas com- 
patíveis com os limites de segurança dos tran- 
sistores. 

Os transistores de germânio conduzem em 
geral a maiores factores de saída do que os de 
silício, por serem menores os valores de Vpr e 
Vcem. Vamos aplicar as equações estabelecidas 
ao dimensionamento duma família de circuitos 
NEM, utilizando o transistor PNP, de germânio, 
ASY32 (Philips). 
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a) Dados 
Transistor 
Lc = — 10 mÃ 
Il; =— 0,2 mÃ 
Vere =— 0,32 V 
Vceu = — 0,2 V 
Icom == — 35 vA (60' €) 


Tensões de alimentação 


Vc =—18BV 
Vccm = — 20 V 
Vcen = — 16 V 
Vos = 18 V 
Vesu = 20 V 
VgBx == 16 V 


Tolerância das resistências (5/0) — y == 0,05 


b) Determinação de Re (Eg. 8) 


R.= 18 ko | 


c) Determinação de Vce (Eq. 9) 


d) Determinação dos parâmetros |, A e B (Egs. 


4, 5 e 6) 
Faz-se VeB=-—-01V e Vcean=0V 
vem | [= 025: mÃ | 
[A=1,051 | 
'B=-0,126 | 


e) Determinação do factor de saída (Eq. 10) 


Fo=15,4 


£) Determinação do maior número de entradas (Ex- 
pressão 11) 


| Mmax == 7 | 
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g) Determinação de Im, Fm, Rm e Rrm (Egs. 3,7,1 e 2) 


Os valores das resistências correspondem aos 
da série de 5º/y mais próximos dos valores cal- 
culados. 


É apresentado como exemplo (Fig. 9) o cir- 
cuito de 3 entradas normalizado. Os factores de 
entrada e de saída inscrevem-se dentro dum 
círculo. 


-18V 


OO 


+18V 


Fig. 9 — Circuito de 3 entradas 


O flip-flop corresponde a 2 blocos de uma 
entrada ligados em circuito fechado. A figura 10 
representa o esquema básico do flip-flop (sem o 
sistema de disparo). 
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ERRATA 


Fig. 10 — Flip-flop 
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— Pág. 161, coluna 1, linhas — 13 e — 12, deve ler-se: 


Vosm (—) —Valor máximo de Vos. | Vosml>|Vegl 


— Pág. 162. Rectificam-se as expressões: 


(4) [= — 
(5) À = | 


Veg 4º Vesn 
Vcem + VeB 
—2D > M+I (1 +29) 


(Ves + VeBn) (1 — 2 7) asc 


(Vc — Voe) (1 — 27) 


— Pág. 162, coluna 1, linha — 10 deve ler-se: 
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1º de Vas 


ERRATA 
ASPECTOS DA ANÁLISE DE CORRELAÇÃO E DE REGRESSÃO 


«Técnica» n.º 359, págs. 481-493 


Em vez de a um nível de probabilidade x 
deve ler-se a um nível de probabilidade 1 — «, nos seguintes locais : 


pág. 
pág. 
pág. 
pág. 
pág. 
pág. 


486, coluna 
486, coluna 
487, coluna 
487, coluna 
489, coluna 
490, coluna 


esg., linha — 13: 
dir., linha— 6; 
esq., linha— 6; 
dir., linha 8; 

dir., linha— 3; 
esq., linha — 12. 
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C. D. DU. 624.31414,5/94 


ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL **) 


I —- Breve nota mensal 


As afluências, em Novembro, ao conjunto de todos 
os aproveitamentos, exceptuados os do Douro, têm 
uma probabilidade de serem excedidas de 39º/y; as 
afluências à totalidade do sistema produtor tém, porém, 
uma probabilidade inferior (28º/,), por efeito do regime 
favorável verificado no Douro. O coeficiente de pro- 
dutibilidade do mês atingiu 1,29, embora as precipita- 
ções tenham sido, como se assinalou, sensivelmente 
inferiores às médias. 


IH — Elementos gerais (GWh) 


a) Mensais Variação 

1965 | 1966 | 9% 
Produção hidráulica (Ph). 3387,1 4690,9, + 38 
Produção térmica (P+)..... 399,3 9,1] — 98 
Produção total (PT). .....| 87794 4700,0 + 24 
Energia recebida de empresas 
não pertencentes ao RNC (Er) 248 31,8 + 3 
Exportações (Ex). . +. .... 1,1) 145)+ 31 
Importações (1) . à mom Rec) AM 3,4 — 99 
Saldo importador (St)... «| 430, q1— 1, 3 — 108 
Consumo em bombagem (Cb) .. 15,8 56 9 + 260 


Produção para con- 

sumos perman. (Pcp)(!).. 
Produção para con- 

sumos não perman. (Pcnp).. 
Total. 
Coeficiente de hidraulicidade 


8758,13946,8| 4 5,0 
460,5, 717,0, + "| 


42844 472907] + 1 
0,48] 1,47 sd 


NOTA 
(1) A produção para consumos permanentes (Pcp) é determinada 


Pep —. Pr + Er + Si ex Cb — Penp 


pela seguinte expressão : 


E —Diagramas de carga dos dias característicos 


4* feira: 


1965 1966 
Produção hidráulica (En) MWh| 16251 17261 
Produção térmica (Pr) MWh 0 0 
Produção total (PT) MWh| 16251 17261 
Trocas com [ Export. (Ex) MWh 0 O 
* Espanha Import. (1) ) Mwh O 0 
Uonsum em bomb. hidroel. (Cb) MWh 690 620 
Prod. para cons. perm.  (Pcp) MWh| 13260 14537 
Prod, para cons, não perm, (Penp) MWh |! 2301 2104 
TOTAL PT + (1-Ex) MW | 15561 16641 
Eotânaiat máx. MW| 850 900 
“o LE otência min, MW 892 405 
2 EE) Utiliz, da ponta horas 18,3 | 18,5 
E E Factor de carga | 0,76 0,77 
PE- Potência máx. MWL To? | Sm. 
io p Potência min. MW “98 319 
e SP Utiliz. da ponta horas 17,5 | 17,8 
Factor de carga 0,73 0.74 
SS es 


NOVEMBRO 


No fim do mês 


Albufeiras: — 


GWh “o (1) 

Alto Rabagão ...... «| 762,1 78,4 
ParIdAiS é sussa val MIS 81,6 
Venda NOVA «sua cus ai 108,8 85,0 
Salamonde . . ... cmol MR 47,5 
Caniçada .co cit. a. : 81,5 95,2 
MEME: ua va sis E Gia 4 5,4 õ,8 
BBB eos é ss E si 181,7 53,6 
Castelo do Bode. +... ,.| 1257 TA 
CM e ea rs E * mu 8,7 44,6 
Lagoa Comprida +. ..... TA (2) 17,5 
Santa Luzia » : wá é é ais 28,2 45,8 
PERORMM oc o wu mg é 5 qa 10,1 78,3 
POMOO ms > srs à dará 4,6(3) 26,0 
Total | com A, Rabagão . . .| 1471,6 69,8 
| sem A, Rabagão . . 708,9 62,4 


NOTAS 


(1) Coeficiente de enchimento em energia das albufeiras. 


(2) Inclui 3.34 GWh armazenados em Vale do Rossim no início de. 
mês e 20 GWh no fim do més, 


(3) Inclui 1,9 GWh armazenados no açude do Poio no início do mês 
e 1,9 GWh no fim do mês, 


(*) Elementos extraídos das estatisticas mensais do Repar- 
tudor Nacional de Cargas (R. N. C.). As produções e 
os consumos das empresas do R. N. C. representam 
cerca de 94º/, dos totais do Pais, 
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Material fabricado nas nossas oficinas, pronto para expedição 


Fabricação de aparelhagem eléctrica de alta e baixa tensão 
Equipamentos electromecánicos completos para centrais e fábricas 
Postos de transformação tipo monobloco 
Quadros blindados, capsulados e em armário 
Armaduras de iluminação de todos os tipos 
Contadores eléctricos e instrumentos de medida 
Telemedida e telecomando 

Instalações elevatórias de água 

Instalações de saneamento e esgotos 
Ferramentas de corte 

Motores Diesel industriais e marítimos 
Fogões eléctricos 

Frigoríficos 

Ventilação e condicionamento de ar 
Aparelhagem de Raios X e electromedicina 
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Durante a germinação, uma semêénte utiliza todos os 
meios ao seu alcance para brotat do solo e poder ver e 
sentir a luz e o calor do Sol. Esta é uma das infinitas 
forças da natureza, que sômente o homem sabe orientar 
de maneira a poder obter o máximo rendimento dela. 

Todos os homens e mulheres, que trabalham para a 
International Harvester, e que se encontram espalhados 
por todo o mundo, ajudam o homem a aproveitar melhor 
essa força da natureza. Através de constantes pesquisas, 
aperfeiçoam a sua técnica, quer no que diz respeito a cada 


INTERNATIONAL 
HARVESTER 


FASSIO, LDA 


uma das partes de uma máquina, quer no que diz respeito 
a produção. Eles produzem tractores e alfaias adequadas 
para cultivar a terra, máquinas para movimentar terras 
e camions para o seu transporte. 

O equipamento 1. H., onde quer que ele seja cons- 
truído, é caracterizado pela sua perfeição, confirmando 
assim o seu lema, de que não são meras ferramentas, 
mas força em si — força aplicável, que satisfaz as neces- 
sidades básicas do homem. 


International Harvester 


põe força nas suas mãos 
20, Rua Jardim do Regedor, 32 


6 LISBOA-2--Telefones 361902-3-4 
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ESTUDO DE UM SISTEMA DE COMANDO AUTOMÁTICO 
PARA O CIRCUITO DE SECAGEM E ABSORÇÃO DE UMA 
FÁBRICA DE ÁCIDO SULFÚRICO 


RESUMO 


Descreve-se um sistema de comando automático desti- 
nado a um processo multidimensional, cuja matriz de 
transferência é não diagonal. O elemento búsico do sistema 
é uma rede correctora que se determina a partir daquela 
matriz e cuja associação ao processo conduz a uma matriz 
de transferência diagonal, É possível, então, formar malhas 
de regulação sem interacção e utilizar controladores con- 
vencionais. Determinan-se, a seguir, os parâmetros destes 
controladores pelo método de Ziegler Nichols. Finalmente, 
estuda-se a estabilidade do sistema pelo metodo do lugar 
geométrico das raizes (diagramas de Evans. 


1 — INTRODUÇÃO 


ENG.º VIRGÍLIO IGREJA 


(Chefe do Dep. «Automação» do 
Centro de Investigação da C. U,F.) 


pelo 


SYNOPSIS 


The non-diagonal transfer matrix of an actual multidi- 
mensional process is determined, From this matriz, a net- 
work qwhich converts that matrix into a diagonal one is 
obtained. The virtually independent loops thus determined 
are then treated by the Ziegler. Nichols"s method in order to 
get the parameters of the P. 1, D. single-loop controllers 
now applicable, At last, the root-locus method is used in 
the stability analysis of the syatem. 


O processo de produção de ácido sulfúrico por contacto é constituído essencialmente pelas 


seguintes fases : 


i) ustulação da pirite; 
ii) purificação dos gases de ustulação ; 
iii) secagem dos gases de ustulação ; 


iv) oxidação catalítica do SO; (que se transforma em SO)); 


v) absorção do SOs. 


As fases de secagem e absorção realizam-se em três torres interligadas, como se indica 
no esquema simplificado da figura 1. Na primeira, chamada torre de secagem, os gases sulfurosos 
são deshumidificados pelo ácido sulfúrico que, finamente pulverizado, atravessa a torre em con- 
tra-corrente. Nas outras duas torres, chamadas torres de absorção, o gás anidrido sulfuroso é com- 


TORRE DE 
ABSORÇÃO 
[1 


Fig. 1 


pletamente absorvido por uma corrente de ácido sulfúrico que, finamente pulverizado, atravessa 
a torre em contra-corrente; desse modo, aumenta-se a concentração do ácido sulfúrico circulante. 
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A necessidade de comando automático nestas fases de fabrico resulta do interesse de se manter 
no tempo, para um maior rendimento do processo, a seguinte situação (*): 


i) concentração aproximadamente constante do ácido à entrada da torre de secagem :95"/01+1º/a; 
ii) concentração aproximadamente constante do ácido à entrada da torre de absorção : 98,5 "/o 
Es 0,1 “o; 
iii) nível de ácido constante (90 + 2) cm na torre de secagem ; 
iv) nível de ácido constante (90 + 2) cm na torre de absorção; 


situação esta que é perturbada pelas seguintes causas aleatórias (entre outras): 


i) variação da quantidade de vapor de água nos gases na torre de secagem (função da tem- 
peratura dos gases à saída dos fornos, temperatura ambiente, etc); 

ii) variação do teor em SO; nos gases na torre de absorção (função do teor em SO: dos gases 
de entrada no convertidor, temperatura de conversão, actividade do catalizador, etc.). 


A anulação dos efeitos destas perturbações sobre a situação desejada, pode fazer-se por mano- 
bras adequadas das válvulas indicadas no esquema da figura 2; contudo, como as variáveis a que 
se refere essa situação estão intrinsecamente ligadas, a manobra de uma qualquer das válvulas 
vai modificar os valores de todas as variáveis definidores da mesma situação, o que torna o comando 
manual do processo muito difícil. 


OH, (vapor) 


| a SOMa (98%) 
“al Sa UA ia a dé A 
ABSORÇÃO 


Pod = 
SO. H; (98%) 


Fig. 2 


O nosso problema consiste em projectar um sistema automático de manobra dessas válvulas. 
Para isso, podemos encarar o processo como um sistema com 4 grandezas de entrada q,, q,, q, e q, 


(caudais determinados pelas válvulas) e 4 grandezas de saída Y,+ Y,» Y;» € Y, (dois níveis e duas 
concentrações). 

O conjunto dos valores das variáveis de entrada pode ser representado por um vector [q,] 
que é transformado por uma matriz P, representativa do processo, num vector de saída [y,] no 
espaço definido pelas variáveis de saída. 


Qy emma O —— Y| [a] 
— da PROCESSO o Cos DJ 
q — — Ya 


Fig. 3 


É possível, então, estabelecer um anel de comando que inclui os instrumentos de medida, os 
órgãos de comando (válvulas) e o controlador segundo o esquema em bloco indicado na fig. 4, onde 
se representam os seguintes vectores: 


[y,] = variáveis de saída ; 


(*) Sem perda de generalidade do método utilizado, consideraremos apenas a torre de absorção I no circuito 
em estudo. 
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[xi] == variáveis de saída traduzidas; 
[xoi] = variáveis de referência traduzidas; 


[4] = [xio — xi] erros; 

[g,] = variáveis de comando do processo ; 
[m;] = variáveis de comando das válvulas; 
[2] = perturbações do processo. 


INSTRUMENTOS 


Fig. 4 


Representando os componentes do anel de comando da fig. 4 por matrizes resulta o diagrama 
da fig. 5. 


A resolução deste problema de comando automático multidimensional, envolve as seguintes 
questões: 


i) Realizabilidade — Com efeito, dado o processo, que não devemos modificar, e dados os 
órgãos de comando (válvulas de comando) e de medida (medidores de nível e de concen- 
tração) cuja liberdade de escolha é muito condicionada, restam-nos livres as característi- 
cas dinâmicas (determinadas pelos valores dos parâmetros livres) do controlador. A pri- 
meira questão consiste em saber se é fisicamente realizável um controlador com tais 
caracte rísticas. 


ii) Estabilidade — Determinada a possibilidade de realização do controlador, é necessário 
determinar a região de estabilidade do sistema, o que vai restringir os intervalos de varia- 
ção dos parâmetros livres do controlador. 


iii) Optimização — Dentro da região de estabilidade e depois de se definir um «índice de 
eficiência» que permita avaliar o comportamento de todo o sistema, há que determinar 
quais os valores dos parâmetros que tornam máximo esse índice, isto é, que tornam 
óptimo, segundo um dado critério, o comando do processo, nas condições impostas ao 


sistema. 
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2 — TEORIA 
2.1 — Esquema geral de controle 


Como se disse atrás, qualquer variável de comando do processo q; influencia todas as variá- 
veis de saída; daqui resulta não adequada a solução que consiste em formar 4 anéis de controle, 
funcionando independentemente, cada um dos quais regulando uma variável y,. Esta solução teria 
a grande vantagem prática de serem imediatamente utilizáveis os controladores existentes no mer- 
cado, da mesma forma que o são as válvulas e os instrumentos de medida necessários ao controle 
do processo. Assim, torna-se necessário realizar um controlador mais complexo que gere os sinais 
adequados [m;] a enviar às válvulas de comando, a partir dos erros [:;], diferenças entre os valores 
fornecidos pelos instrumentos de medida e os valores desejados para as variáveis de saída. Contudo, 
com o fim de beneficiarmos daquela enorme vantagem de ordem prática, vamos basear o nosso 
estudo num método que consiste em neutralizar a influência recíproca existente entre as variáveis de 
entrada e de saída do sistema, pela inclusão duma rede correctora na cadeia activa do anel de 


comando, como se indica na figura 6. 
[2 ) 
PROCESSO 


CONTROL ADORE 5 
CONVENCIONAIS 


INSTRUMENTOS 


Fig. 6 


Nestas condições, cada um dos controladores convencionais vai actuar apenas em uma variável 
de saída y, embora, para isso, tenha em geral que actuar em todas as variáveis de comando q, . 
De qualquer modo, cada controlador irá comportar-se como se a ele estivesse associado apenas um 
anel de controle, uma vez que a partir dum sinal de erro s; gerará um sinal que actuará únicamente 
a variável de saída y,. Daqui resulta ainda a possibilidade de utilização dos métodos de análise 
dos sistemas unidimensionais, nomeadamente, os de Evans e Ziegler-Nichols para o estudo da esta- 
bilidade de cada um dos anéis. 

Assim, o esquema geral de controle que vamos utilizar é o indicado na figura 7. Ele servirá 
de base à dedução das equações teóricas de comando que faremos a seguir. 


do mm 
E a 
E PROCESSO 
ê 
Z, 
V 
1 A E, 5 my q, | | y, 
Gt. SS) - =| válvula Tc S) 
RED | y 
toa AE Com TJ ; [váluto 7] 0x4 
É s CORREC- PROCESSO 
e dE la 20 Srs 
“oe dO E Cont a js - - q — 
= [nstrum 1] — 
[instrum 3] 
Fig. 7 
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2.2 — Equações básicas 


O modelo físico descrito pelo esquema geral de controle estabelecido anteriormente supõe 
os componente lineares e invariantes no tempo, isto é, supõe que o comportamento dinâmico dos 
mesmos pode ser descrito por um sistema de equações diferenciais lineares de coeficientes constantes. 

Utilizando a transformada do Laplace para a determinação das funções de transferência cor- 
respondentes aos componentes indicado na figura 7 e, por outro lado, utilizando a linguagem matri- 
cial para a dedução das equações básicas de funcionamento, obtém-se o seguinte esquema a partir 
do anterior: 


Fig. 8 


Neste esquema C, V e I representam matrizes diagonais de transferência correspondentes 
aos conjuntos dos controladores, das válvulas e dos instrumentos de medida; P e H representam 
as matrizes correspondentes ao processo que transformam respectivamente os vectores das variáveis 
de comando [q,] e das perturbações [z,] nas parcelas que constituem o vector de saída [y,]; [x] é o 
vector de saída traduzido em termos susceptíveis de comparação com o vector de entrada [x,;] de 
modo a fornecer o vector erro [:]; finalmente, [s;] e [m,] são, respectivamente, os vectores de en- 
trada e saída da rede correctora descrita pela matriz R. Os elementos de todas estas matrizes 
e vectores são funções da variável complexa s == 2 + jw em consequência da aplicação da transfor- 
mada de Laplace ao sistema de equações diferenciais que descrevem o funcionamento do sistema. 

A partir da figura 8, obtêm-se as seguintes equações matriciais : 


ly) =P [lg] + Hz] (1) 
[q] = VRC [5] (2) 
[5] = |x] — [x] (3) 
[x] = [yi] (4) 

donde resulta, fazendo G = PVRC: 
ly]=Gl[]+HIz] (5) 

e, se a matriz (E + Gl) —! existir: 

Ly] = (E + GD”! G[x] + (E + GD”! H [2] (6) 


sendo E a matriz unidade. 


Destas expressões é possível determinar, por ensaios, V, |, Pe H; da matriz C, diagonal, 
conhecemos os limites e intervalos de variação dos seus elementos (*); resta-nos determinar a 
matriz R de acordo com as condições de realizabilidade, estabilidade e optimização que se indicam 
a seguir. 


(*) Utilizaremos controladores convencionais P. 1. D, (controle proporcional, integral e derivativo). 
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2.3 — Realizabilidade 


Para que o sistema seja realizável, tal como o concebemos, devem ser satisfeitas as seguintes 
condições : 


i) existência da matriz inversa (E + GI)”!, em virtude da equação (6); 
ii) carácter diagonal da matriz G == PVRC, em virtude da equação (5) e para que não haja 
influência de +, sobre y,, para k & i; 
iii) valores teóricos razoáveis ligados aos componentes e variáveis associados à rede correctora. 


2.4 — Estabilidade 
Da equação (6) conclue-se que a equação característica do nosso sistema é: 


IE+GI|=0 (7) 


onde G ==PVRC. Como se sabe, todas as raízes desta equação em s==7 + jw devem possuir partes 
reais negativas para que o sistema seja estável. Esta condição restringe os valores dos elementos 
&; de Ger, de R. Por outro lado, em virtude do carácter diagonal da matriz global de transfe- 
rência tudo se passa como se o nosso sistema fosse constituído por 4 anéis de controle independen- 
tes; isso permite-nos utilizar, para cada um deles, o método do lugar geométrico das raízes (método 
de Evans) na escolha duma boa região de estabilidade. Desta condição resultam novas limitações 
para os valores de g;, de r; e dec, de €. 


2.5 — Optimização 


As condições anteriores de realizabilidade e estabilidade não conduzem a um conjunto único 
de valores dos parâmetros livres do sistema. Com efeito, permitem-nos escolher, apenas, no espaço 
dos parâmetros livres, uma região permissível 2, para os valores desses parâmetros que garantem 
estabilidade ao sistema. Existe porém uma região “',; C 4; para a qual o sistema tem um compor- 
tamento que consideramos «óptimo» em face do tipo de perturbações a que está sujeito. Um dos 
critérios utilizados na avaliação do comportamento do sistema é o do erro quadrático médio, através 
do método de Phillips. Não faremos o respectivo estudo neste trabalho. 


2.6 — Aplicação ao problema em estudo 


Do esquema da fig. 8, podemos determinar experimentalmente P, V, I e as características de [z;]. 

Por outro lado, conhecemos os intervalos de variação dos parâmetros livres dos controladores 
que definem C e o vector de entrada [xoi]. O nosso problema consiste na determinação de R e dos 
parâmetros de C, a partir apenas destes dados, de tal modo que o sistema seja realizável e estável 
nos sentidos dados atrás. 

Na resolução deste problema seguiremos uv método seguinte: 


i) Determinação experimental das funções de transferência parciais que constituem os ele- 
mentos das matrizes P, Ve 1. 


ii) Determinação dos elementos r,, de R, de modo que G' =IPVR seja diagonal e g', sejam 


—Lg s 
Kg e 


Cgs+1 
serem nulas algumas das constantes Ku, Lg e Cs), ou dele se possam aproximar, com as 
constantes relacionadas, através das fórmulas de Cohen-Coon, com os parâmetros permis- 
síveis para Xp, T+ e Ti, dos controladores convencionais. 

iii) Verificação das condições de realizabilidade i), ii), iii) indicadas em 2.3. 

iv) Verificação da estabilidade dos anéis de controle pelos respectivos diagramas de Evans. 


do tipo g",;(s) = (tempo morto e constante de tempo, podendo eventualmente, 
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3 — PROJECTO DO SISTEMA CONTROLADOR 


3.1 — Determinação experimental da matriz de transferência do processo 


3.1.1 — Esquema em bloco do processo 


A fim de prepararmos os ensaios necessários à determinação da matriz P, consideremos o 
esquema da figura 9, onde se assinalam todos os elementos que interessam àqueles ensaios. Neste 
esquema tem-se: 


9,2% A22 33 


Fig. 9 


v,— válvula de admissão de água; 


Tm válvula de passagem de ácido a 98º/ do circuito de absorção para o circuito de 
secagem ; 


v— válvula de passagem de ácido a 95/o do circuito de secagem para o circuito de 
absorção ; 


v, — válvula de expedição de ácido a 98 0; 
q; — caudal em volume de água comandado pela válvula v, (= 0,25 m'h-!); 
qv, — caudal em volume de ácido a 98 “/o comandado pela válvula v, (=7 m'h-!); 
qv, — caudal em volume de ácido a 95 */y comandado pela válvula vs (=7 m'h-!); 


qv, — caudal em volume de ácido a 98º/o de expedição, comandado pela válvula 
v; (=0,6 m*h”!); 


qv — caudal em volume de água absorvida na torre de secagem (= 0,25 m'h !); 

q, — caudal em massa de anidrido sulfúrico absorvido na torre de absorção; 
am caudal em volume do ácido de saída da torre de secagem com a concentração y,; 
q, Y, — caudal em volume do ácido de entrada na torre de secagem, com a concentração Y,? 
Gy, caudal em volume de ácido de entrada na torre de absorção com a concentração Y,5 
Hd = caudal em volume do ácido de saída da torre de absorção com a concentração y, ; 

qs, — caudal em volume de ácido da bomba B, (- 30 m'h-!); 

Ja, — caudal em volume de ácido da bomba B, (=30 m*h”!); 

Rf, — refrigerante do circuito de secagem (conduta com cerca de 260 m); 

Rt, — refrigerante do circuito de absorção (conduta com cerca de 130 m); 


Rf, — refrigerante do circuito de expedição (conduta com cerca de 10 m); 
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y,— nível de ácido na torre de secagem (= 90 cm); 
Y, — concentração do ácido à entrada da torre de secagem (= 95º) 
Y, — nível do ácido na torre de absorção (= 90 cm); 


Y, — concentração do ácido à entrada da torre de absorção (98,5 /9). 


Da figura 9 resulta o seguinte esquema em bloco geral do processo (fig. 10). 
3.1.2 — Dedução teórica do tipo de matriz de transferência 


Como base no esquema da figura 10 vamos construir um modelo físico que descreva o pro- 
cesso com uma aproximação razoável nas condições estacionárias pretendidas. Para isso, note-se, em 
primeiro lugar, que as grandes perdas de carga se localizam, quer no circuito de secagem quer no 
circuito de absorção, nos refrigerantes Rfi e Rfs, já que estes constituem cerca de 90º 0 do compri- 
mento total das condutas. Por outro lado e como se indicou já na secção anterior, os caudais gB1 e 
gB: são da ordem dos 30 m'h-! e os caudais qvz e qvs da ordem dos 7 m'h”-!. Considerando 
variações máximas de + 15º/o para os caudais de comando qvs e qvs, elas representam, em relação 
a qBi e qB>, sômente + 3º/y dos valores destes caudais; daqui resultam, em relação às perdas de 
carga em Rfy e Rfs, variações desprezáveis. 

Se acrescentarmos que as velocidades das bombas B; e Bs são constantes, podemos concluir 
ser válida a hipótese que faremos, de serem qBi e qBz constantes nas condições normais de regu- 
lação. Outra hipótese de partida é a seguinte: nos nós 1 e 3, a mistura faz-se perfeitamente e sem 
contracção de volume (!). Por outro lado, vamos admitir que as variações de densidade depois das 
misturas, nos nós 1 e 3, são desprezáveis em relação às correspondentes variações de concen- 
tração (*). Finalmente, supomos o fluido incompressível e as condutas indeformáveis. 

Partindo destas hipóteses vamos agora deduzir o sistema de equações diferenciais que dará 
origem à matriz de transferência P definida pela equação matricial (1) supondo apenas variações das 
variáveis de comando q, : 


Ly; (s)] = P(s). [g, (s)] (8) 


Da figura 10 obtém-se o esquema simplificado da figura 11 que interessa ao nível y, da torre 
de secagem. 


Fig. 11 


Pelas hipóteses feitas, resulta para os nós 1 e 3: 


9B, = qu + qvua= qv; + qv 
donde 
qu =gqB,— qvs 


a = q9B, — qva 


(!) Em misturas de ácido a 100º/, com água, a 15ºC, o factor de contracção é 0,957. 

(2) Dos valores indicados na literatura obtivemos, para 15ºC e concentrações próximas de 95º/q e 98º/9, resul- 
tados inferiores, respectivamente, a 10º/y e 5º/o, para os valores máximos das relações entre as taxas de variação da 
densidade e da concentração. 
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Substituindo estes valores na equação de balanço volumétrico 


dy, = 
dt 


A, qv + Ga + qm — 


onde A, é a secção da torre cilíndrica de secagem, vem: 
d 
A, “a = qv + qv, + qu — qv; (9) 


Se ya, QV, Vi, QVa € qv; representarem desvios em relação a valores de referência y,9, QVo, 
dVio + Voo + Va » Obtém-se de (9) por aplicação da transformada de Laplace : 


A,sY.(s) = Qv (s) + Qv, (5) + Qvs(s) — Qv; (5) 


donde resultam as funções de transferência seguintes, elementos da matriz procurada P: 


P. (5) Ys (s) 1 
s) = - E quere 
= Qv, (s) A,s 
Ya (s) 1 
P E mm E eme 
ne (8) Qvo (s) A,s 
(10) 
Ya (s) 1 
P = í— 0 = — — 
e 49) Qva (s) A,s 
Ya (s) 
P + qu 
24 (s) Qu, (8) 


Anãlogamente, em relação ao nível y, temos: 


dB, = q, Te av, = av, av, Uy 


Fig. 12 


donde 


Ta = 96, — av, — qv, 
Substituindo estes valores na equação de balanço volumétrico 


dy, 
2? dt Las — Ig 
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onde As é a secção da torre cilíndrica de absorção, vem: 
dy, 


A, a 4 Va avi (11) 


de y, € q, representarem desvios em relação aos valores de referência Yao Ea é 


obtém-se de (11) pela aplicação da transformada de Laplace: 
A, sY, (s) = Qv, (s) — Qv, (s) — Qv, (s) 


donde resultam as funções de transferência seguintes, elementos da matriz procurada : 


q Re 26 nt 

É cdi E Eae E 

Lada a = e 
P,() = Es 


Consideremos agora o esquema simplificado seguinte, derivado do da fig. 10, para a dedução 
das equações que interessem à concentração y,: Nesta figura P, e T, designam, respectivamente, O 
processo e a torre de secagem, y, a concentração que resulta para o ácido de entrada depois da 
secagem, y, a concentração do ácido q di da torre de absorção, sendo conhecidos já os outros símbolos. 


Por definição de concentração, tem-se, 


ED d, Y 
qa d, + IJ Ea 4, 


onde d, designa a densidade do ácido à concentração y, 


y; == 


Fazendo o balanço mássico de ácido na torre, resulta a equação: 


d 
Veda e = (Qua + qu) dyys— Qu da Yu 


donde resulta, pela substituição do valor y.: 


d (qr + qo) quad do 
Vo dy a + QudyuYu = Sue Se A. d di Sa 
t Jo ] + dy A qu 
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Substituindo nesta expressão os valores de q,, e qu dados por: 


Jg — dp, — Qua 


Ju = qo, — Ty 


obtém-se 


d 
Vody e + GmduYu— Ig du = É (dy Goyr Togr Ya) (13) 


com 
(qui — Iv; + Gy) (Gr, — dog) ddçyi 


(q, — Ig) da ++; a 


É (qo Joy Toy Y4) = 


Faremos algumas aproximações em (13). Em primeiro lugar desprezaremos a pequena parcela 
resultante do desvio de dy, NO 1.º membro, em relação à grande parcela (gp, + dog) du yum. Por 


outro lado e supondo pequenos desvios de g,, dy” Ivy € Y1, VAMOS linearizar o 2.º membro de (13) 


considerando apenas as primeiras parcelas do respectivo desenvolvimento em série de Taylor. 
Obtemos, assim, uma equação diferencial do tipo: 


V. dy 
dt 


+ IgJu=C +Cq, + Co t+ Cos + Coys (14º) 


onde C, €,, €,, €,, e €, são constantes. 
. — A “ 
Supondo Yy,,, ds Jo € Ty desvios em relação a valores de referência Yuo! Too Togo! 


Togo * Yi” obtém-se de (14º) pela aplicação da transformada de Laplace: 


1 
VD eV+ra)=DCT+ECOQ (O) +C, O (8) + Co Qu; (8) + CY, (8) (15) 
Em relação ao nó N, da figura 13, tem-se: 


IT d,, Yu + Jyg d, Y3 


y, = (16) 
Ja d, d, d, 
Linearizando esta expressão, tal como fizemos em (14), vem: 
q = E Yu É CR ERR +- CY, (17) 
cuja equação transformada é: 
À au 
M (5) = e Cr + ni po (s) + Cy, Qua (s) + Co E (18) 


Retomando a figura 13, consideremos o refrigerante R, - À variável de saída é a mesma de 
entrada: a concentração Y,; como a conduta é longa, porém, há um atrazo sensível t, que dá ori- 
gem à resposta y, (t — t,) a uma entrada Y, (t). Representando a concentração Y,* no nó N |" por 
Y, Ng e no nó N , por +, , podemos escrever: 


YO = yn, €—t) (19) 


equação cuja transformada é; 
E (8) = et “YAN, (s) (20) 


178 TECNICA N. 364 


Estamos agora em condições de determinar a equação diferencial transformada global corres- 
pondente à cadeia de elementos em série indicada na figura 14. Assim, a partir de (15), (18) e (20) 
obtém-se : 


C, Cy erts É es a 18 * Co Cu —t s Co Cu —ty 8 
“0 (= ra ) Te CS ra E a + 
—t; 8 
Si Q,, (8) + A et Q, (8) + CH SQ, (9) + Cetisy,(s) (21) 
sV. + Ie, E sV, + By ' i 


Supondo em (21) apenas variações de Q,, (s), Q,, (s), Q,, (s) e Y, (s), podemos escrever: 


C, C Co, Cm Co, Cu 
Y69 (1 ee ) = et Q (D+ LL ensQ, (9) + 
Ê sV, + de, sV, + dB, 4 sV.+ dy o 


+ Ca Q, (8) + Ce Y, (6) (22) 


Consideremos, agora, o esquema da figura 14 (extraído da figura 10) a fim de deduzirmos 
as equações que interessam à variável concentração Y3. 


Os símbolos P, e T, designam, respectivamente, o processo e a torre de absorção; y, repre- 
senta a concentração resultante do processo de absorção. Seguindo um raciocínio semelhante ao uti- 
lizado em relação a y, obtém-se para y, a seguinte equação: 
Ea € , aa; Cu Cy € 
Va + ds Va + q, Va + q 


XY, (9) (1 — Q,, (5) + 


+ Cu emtis Q,, (s) + C,ertrs Y, (s) (23) 


Obtidas as equações (22) e (23), podemos finalmente formar o sistema (24) do qual resultam 
as expressões de Y, (s) e Y, (s) em função de Q, (s), Q, (s), Q,, (s) e Q, R (s): 


——— 


C, Cy erts Co Cu erts Cu Ce 
hi (1 ú sV, + qu Jr Va Qu (s) + Vs + ds, Qus (s) + Cy q Qu, (s) + 
+ Cert4sY, (s) 
Ys (s) (1 RE — Ca et Q,. (5) + a et Q,, (s) + 
SVa É, / Va + du, , SVa + du, ) 
+ Co?" Qu (8) + C et Y, (8) (24) 
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Resolvendo estas equações obtêm-se as seguintes expressões gerais para Yi(s)eY; (s): 


, Ei tis EO, | 
| (s) — Vas, e E 
e 


1 


À — — 
| sV, + q, 


— tis — tis " ca ; vo" 3 
e + 9, e Q,, (s) + Coe — aaa: 


6 Cu — tas 
es ii sSV +q 
| E Va + Im ) ( á so 


ti C,C'v; nr ot + to)s | 
“ia ” CC o —W 
TE SPSS. 


— (tit-ty)s 
— e 


4 A (s) ap 
E sd Es e tis ; 
Vs a IB, ) 


ta E mA E: 70 e ltitt)s Q,, (8) EE 
a 2 C, Cu 
sV, + dp, 


—L(tHt)s 


Qu, (8) + 


(25) 
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